Improvements to the APBS biomolecular solvation software suite

The Adaptive Poisson-Boltzmann Solver (APBS) software was developed to solve the equations of continuum electrostatics for large biomolecular assemblages that has provided impact in the study of a broad range of chemical, biological, and biomedical applications. APBS addresses three key technology challenges for understanding solvation and electrostatics in biomedical applications: accurate and efficient models for biomolecular solvation and electrostatics, robust and scalable software for applying those theories to biomolecular systems, and mechanisms for sharing and analyzing biomolecular electrostatics data in the scientific community. To address new research applications and advancing computational capabilities, we have continually updated APBS and its suite of accompanying software since its release in 2001. In this manuscript, we discuss the models and capabilities that have recently been implemented within the APBS software package including: a Poisson-Boltzmann analytical and a semi-analytical solver, an optimized boundary element solver, a geometry-based geometric flow solvation model, a graph theory based algorithm for determining p$K_a$ values, and an improved web-based visualization tool for viewing electrostatics

Automated Detection and Segmentation of Synaptic Contacts in Nearly Isotropic Serial Electron Microscopy Images

We describe a protocol for fully automated detection and segmentation of asymmetric, presumed excitatory, synapses in serial electron microscopy images of the adult mammalian cerebral cortex, taken with the focused ion beam, scanning electron microscope (FIB/SEM). The procedure is based on interactive machine learning and only requires a few labeled synapses for training. The statistical learning is performed on geometrical features of 3D neighborhoods of each voxel and can fully exploit the high z-resolution of the data. On a quantitative validation dataset of 111 synapses in 409 images of 1948×1342 pixels with manual annotations by three independent experts the error rate of the algorithm was found to be comparable to that of the experts (0.92 recall at 0.89 precision). Our software offers a convenient interface for labeling the training data and the possibility to visualize and proofread the results in 3D. The source code, the test dataset and the ground truth annotation are freely available on the website http://www.ilastik.org/synapse-detection

Semi-Automated Reconstruction of Neural Processes from Large Numbers of Fluorescence Images

We introduce a method for large scale reconstruction of complex bundles of neural processes from fluorescent image stacks. We imaged yellow fluorescent protein labeled axons that innervated a whole muscle, as well as dendrites in cerebral cortex, in transgenic mice, at the diffraction limit with a confocal microscope. Each image stack was digitally re-sampled along an orientation such that the majority of axons appeared in cross-section. A region growing algorithm was implemented in the open-source Reconstruct software and applied to the semi-automatic tracing of individual axons in three dimensions. The progression of region growing is constrained by user-specified criteria based on pixel values and object sizes, and the user has full control over the segmentation process. A full montage of reconstructed axons was assembled from the ∼200 individually reconstructed stacks. Average reconstruction speed is ∼0.5 mm per hour. We found an error rate in the automatic tracing mode of ∼1 error per 250 um of axonal length. We demonstrated the capacity of the program by reconstructing the connectome of motor axons in a small mouse muscle

Method: automatic segmentation of mitochondria utilizing patch classification, contour pair classification, and automatically seeded level sets

<p>Abstract</p> <p>Background</p> <p>While progress has been made to develop automatic segmentation techniques for mitochondria, there remains a need for more accurate and robust techniques to delineate mitochondria in serial blockface scanning electron microscopic data. Previously developed texture based methods are limited for solving this problem because texture alone is often not sufficient to identify mitochondria. This paper presents a new three-step method, the Cytoseg process, for automated segmentation of mitochondria contained in 3D electron microscopic volumes generated through serial block face scanning electron microscopic imaging. The method consists of three steps. The first is a random forest patch classification step operating directly on 2D image patches. The second step consists of contour-pair classification. At the final step, we introduce a method to automatically seed a level set operation with output from previous steps.</p> <p>Results</p> <p>We report accuracy of the Cytoseg process on three types of tissue and compare it to a previous method based on Radon-Like Features. At step 1, we show that the patch classifier identifies mitochondria texture but creates many false positive pixels. At step 2, our contour processing step produces contours and then filters them with a second classification step, helping to improve overall accuracy. We show that our final level set operation, which is automatically seeded with output from previous steps, helps to smooth the results. Overall, our results show that use of contour pair classification and level set operations improve segmentation accuracy beyond patch classification alone. We show that the Cytoseg process performs well compared to another modern technique based on Radon-Like Features.</p> <p>Conclusions</p> <p>We demonstrated that texture based methods for mitochondria segmentation can be enhanced with multiple steps that form an image processing pipeline. While we used a random-forest based patch classifier to recognize texture, it would be possible to replace this with other texture identifiers, and we plan to explore this in future work.</p

Imaging Transient Blood Vessel Fusion Events in Zebrafish by Correlative Volume Electron Microscopy

The study of biological processes has become increasingly reliant on obtaining high-resolution spatial and temporal data through imaging techniques. As researchers demand molecular resolution of cellular events in the context of whole organisms, correlation of non-invasive live-organism imaging with electron microscopy in complex three-dimensional samples becomes critical. The developing blood vessels of vertebrates form a highly complex network which cannot be imaged at high resolution using traditional methods. Here we show that the point of fusion between growing blood vessels of transgenic zebrafish, identified in live confocal microscopy, can subsequently be traced through the structure of the organism using Focused Ion Beam/Scanning Electron Microscopy (FIB/SEM) and Serial Block Face/Scanning Electron Microscopy (SBF/SEM). The resulting data give unprecedented microanatomical detail of the zebrafish and, for the first time, allow visualization of the ultrastructure of a time-limited biological event within the context of a whole organism

A Computational Framework for Ultrastructural Mapping of Neural Circuitry

Circuitry mapping of metazoan neural systems is difficult because canonical neural regions (regions containing one or more copies of all components) are large, regional borders are uncertain, neuronal diversity is high, and potential network topologies so numerous that only anatomical ground truth can resolve them. Complete mapping of a specific network requires synaptic resolution, canonical region coverage, and robust neuronal classification. Though transmission electron microscopy (TEM) remains the optimal tool for network mapping, the process of building large serial section TEM (ssTEM) image volumes is rendered difficult by the need to precisely mosaic distorted image tiles and register distorted mosaics. Moreover, most molecular neuronal class markers are poorly compatible with optimal TEM imaging. Our objective was to build a complete framework for ultrastructural circuitry mapping. This framework combines strong TEM-compliant small molecule profiling with automated image tile mosaicking, automated slice-to-slice image registration, and gigabyte-scale image browsing for volume annotation. Specifically we show how ultrathin molecular profiling datasets and their resultant classification maps can be embedded into ssTEM datasets and how scripted acquisition tools (SerialEM), mosaicking and registration (ir-tools), and large slice viewers (MosaicBuilder, Viking) can be used to manage terabyte-scale volumes. These methods enable large-scale connectivity analyses of new and legacy data. In well-posed tasks (e.g., complete network mapping in retina), terabyte-scale image volumes that previously would require decades of assembly can now be completed in months. Perhaps more importantly, the fusion of molecular profiling, image acquisition by SerialEM, ir-tools volume assembly, and data viewers/annotators also allow ssTEM to be used as a prospective tool for discovery in nonneural systems and a practical screening methodology for neurogenetics. Finally, this framework provides a mechanism for parallelization of ssTEM imaging, volume assembly, and data analysis across an international user base, enhancing the productivity of a large cohort of electron microscopists

A Genetically Encoded Tag for Correlated Light and Electron Microscopy of Intact Cells, Tissues, and Organisms

Electron microscopy (EM) achieves the highest spatial resolution in protein localization, but specific protein EM labeling has lacked generally applicable genetically encoded tags for in situ visualization in cells and tissues. Here we introduce “miniSOG” (for mini Singlet Oxygen Generator), a fluorescent flavoprotein engineered from Arabidopsis phototropin 2. MiniSOG contains 106 amino acids, less than half the size of Green Fluorescent Protein. Illumination of miniSOG generates sufficient singlet oxygen to locally catalyze the polymerization of diaminobenzidine into an osmiophilic reaction product resolvable by EM. MiniSOG fusions to many well-characterized proteins localize correctly in mammalian cells, intact nematodes, and rodents, enabling correlated fluorescence and EM from large volumes of tissue after strong aldehyde fixation, without the need for exogenous ligands, probes, or destructive permeabilizing detergents. MiniSOG permits high quality ultrastructural preservation and 3-dimensional protein localization via electron tomography or serial section block face scanning electron microscopy. EM shows that miniSOG-tagged SynCAM1 is presynaptic in cultured cortical neurons, whereas miniSOG-tagged SynCAM2 is postsynaptic in culture and in intact mice. Thus SynCAM1 and SynCAM2 could be heterophilic partners. MiniSOG may do for EM what Green Fluorescent Protein did for fluorescence microscopy